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' Abstract. 

We present FCRAO maps as well as combined BIMA and FCRAO maps of the high density molecular emission towards the 
■ CS core in the L673 region. With the FCRAO telescope we mapped the emission in the CS (2-> 1), C^**S (2-» 1), HCO+ 

' (1— » 0), and H'^CO^ (1— > 0) lines. The high-density molecular emission, which arises from a filamentary structure oriented 
in the NW-SE direction, shows clear morphological differences for each molecule. We find that HCO^ has an extremely high 
04 ' optical depth, and that the H'^CO^ emission is well correlated with submm sources. The BIMA and FCRAO combined maps 

^ I recover emission from structure previously undetected or marginally detected, and show an overall aspect of a filamentary 

structure connecting several intense clumps. We found a total 15 clumps in our combined data cube, all of them resolved at our 
angular resolution, with diameters in the 0.03-0.09 pc range. Their estimated masses range between 0.02 and 0.2 Mq, except 
for the largest clump, which has a mass of ~ 1.2 Mq. We find a clear segregation between the northern and southern region 
of the map: the northern section shows the less chemically evolved gas and less massive but more numerous clumps, while the 
' southern region is dominated by the largest and most massive clump, and contains the more evolved gas, as traced by emission 

, of late-time molecules. We find that the derived clump masses are below the virial mass, and that the clump masses become 

closer to the virial mass when they get bigger and more massive. This supports the idea that these clumps must be transient, and 
that only the more massive ones are able to condense into stars. The clumps we detect are probably in an earlier evolutionary 
stage than the "starless cores" reported recently in the literature. Only the most massive one has properties similar to a "starless 
Q ' core". 

C/3 . Key words. ISM: individual objects: L673 — ISM: abundances — ISM: clouds — ISM: molecules — Radio lines: ISM 

^ ' Stars: formation 
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^ 1 . Introduction lar the systematic difference in size and separation between the 

' emission peaks of both molecules. 
Dense cores of molecular clouds have been extensively used 

to study the star formation process. The emission of several To explain these results, Taylor et al. CI 996" 1998"! devel- 

molecules known to be good tracers of high density molecu- oped a chemical model in which high density condensations, 

lar gas has been used to study their properties, such as mass or dense cores, are formed by clumps of different masses, 

and density, temperature, and velocity distributions. However, ages, sizes, and densities but with sizes of $ 0.1 pc. These 

it was apparent from the first surveys (Zhou et al. .l989. Myers clumps would be unresolved in single-dish angular resolution 

et al. I1991> that there existed large discrepancies between the observations such as the ones listed in the previous paragraph, 

emission of some of these molecules, such as CS, NH3 or Most of the clumps would disperse before NH3 abundances 

HCO^, which needed to be clarified in order to obtain the ac- build up to significant levels, though these clumps would con- 

tual distribution of the high density gas in these sources. Pastor tain substantial CS, so its emission should be observable. A 

et al. (TWT) and Morata et al. (19971 earned out a system- few clumps, those sufficiently long-lived, or in a more ad- 

atic comparison between the emission of the CS (7=1— >0) and vanced stage of physical and chemical evolution because of be- 

NH3 (J, ^r)=(l,l) transitions, under similar conditions of angu- ing denser or more massive, would form a significant amount 

lar resolution, which confirmed these discrepancies, in particu- of NH3, while the CS abundance decreases with time. These 

clumps would possibly continue their evolution to eventually 
form stars. Moreover, Taylor et al. (119961 11998> found that 
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"late-time " molecules (showing more compact emission like 
NH3), according to the time at which these molecules reached 
their peak fractional abundance. 

In order to test the Taylor et al. chemical model, and in 
particular to study the small-size structure of the dense cores, 
we carried out high angular resolution observations with the 
BIMA telescope towards the L673 molecular cloud of sev- 
eral transitions corresponding to molecules of both families of 
species (Morata, Girart & Estalella 2003', hereafter MGE03). 
The region selected had been previously mapped at an angu- 
lar resolution of ~ 1' in CS (2^ 1) (Morata et al. lT997l and 
NH3 (J,K)-{1,1) (Sepiilveda et al. in preparation). It has no 
signs of star formation and is located ~ 9'5 to the southeast of 
IRAS 19180-Hlll6.The distance to the L67 3 cloud is ~ 300 pc 
(Herbig & Jones IT^ 

MGE03 detect emission in the CS (7=2^1), N2H+ 
(7=1^0), HCO+ (7=1^0) and SO (Jk = 23 ^ I2) lines. 
The high angular resolution interferometric observations reveal 
a much clumpier medium than the lower resolution single-dish 
observations. Several clumps of size $ 0.08 pc are detected 
distributed all over the region encompassed by the single-dish 
observations. However, the BIMA interferometer only detects 
9-12% of the CS (2->l) emission. ModelHng the filtering ef- 
fect of BIMA, MGE03 show that a clumpy and heterogeneous 
medium could explain this effect. 

MGE03 find that the BIMA observations can be fitted by 
a model with a density at which collapse is halted of «/ = 
5 X 10^ cm"-^. A tentative classification of the studied clumps 
according to the stage of chemical evolution indicated by the 
molecular abundances is proposed. However, the model cannot 
explain the observed low HCO^ abundances with respect to the 
predicted values obtained at two positions, but this is possibly 
due to strong self-absorption in the HCO^ (1— » 0) line. 

The results of MGE03 show that although the interferomet- 
ric observations lose information on the total emission of the 
region, they are in agreement with the results obtained with 
single-dish telescopes, suggesting that interferometric obser- 
vations are useful for studying with high angular resolution 
the chemical evolution of the inner structures of starless cores. 
However, in order to avoid the filtering effect of interferome- 
ters, array observations should be combined with single dish 
observations. This would give high resolution maps of the full 
emission. 

In this paper we present the study of the spatial distribu- 
tion of the molecular emission in L673 from the maps obtained 
with the 14-m FCRAO telescope, and from the maps obtained 
by combining the BIMA data of MGE03 and the FCRAO maps 
presented in this paper. In Sect. 2 we describe the FCRAO 
observations and the BIMA and FCRAO combination proce- 
dure. In Sect. 3 we describe the morphological properties of 
the molecular emission. In Sect. 4 we analyze the properties of 
the molecular cloud observed. 



Table 1. Lines Observed with the 14-m FCRAO telescope 







V 


FWHM 


Av 




Molecule 


Transition 


(GHz) 


Beam 


(km/s) 


(K) 


H"CO+ 


(1^0) 


86.75434 


59" 


0.17 


0.04 


HCO+ 


(1-^0) 


89.18852 


58" 


0.16 


0.11 




(2^ 1) 


9641295 


53" 


0.15 


0.07 


CS 


(2^ 1) 


97.98097 


52" 


0.15 


0.14 



2. Observations 

2.1. FCRAO Observations 

The observations were carried out in 2003 January with the 
14-m Five College Radio Astronomy Observatory (FCRAO) 
telescope using the SEQUOIA (Second Quabbin Observatory 
Imaging Array) 16 beam array receiver We used the dual chan- 
nel correlator (DCC), which provided a bandwidth of 25 MHz 
and a spectral resolution of 24.4 kHz. The DCC was config- 
ured to observe the 4 molecular transitions using 3 different fre- 
quency setups: one for HCO^, one for CS, and one for H'^CO^ 
and C^'^S, which were observed simultaneously. The main- 
beam efficiency at the frequencies observed is 7]^^ ^ 0.50. 
At 89 and 98 GHz the half-power beam width is 58" and 
52", respectively. The weather conditions during the obser- 
vations were excellent, providing averaged system tempera- 
tures of ~120 K. The observations were performed using the 
On-The-Fly (OTF) mapping technique, covering a region of 
8' X 8' for HCO+ and CS and of 6' x 6' for H'^CO^ and 
C^'^S. The center of the maps was Qf(72000) = 19^20'"5Z20; 
5(72000) = 11°14'45"5. The selected reference position was 
30' east of the center position. A summary of the observed 
molecular line transitions is given in Table ^ The OTF maps 
were regridded to have a cell size of 20". The line intensities 
are given in main-beam brightness temperature, T^b = 7'^/'7mb- 

2.2. Combination procedure for the FCRAO and BIMA 
data 

The BIMA interferometer detects only 9-12% of the CS 
(2— > 1) emission. Therefore, in order to recover the missing 
flux of the interferometer data we combined the CS (2— > 1) 
and HCO+ (1-^0) data obtained with the FCRAO presented in 
this paper with the BIMA data from MGE03. The first step 
to do this is to convert the FCRAO antenna temperature to 
a flux scale. We used a 43 Jy K"' conversion factor for the 
two transitions (M. Heyer, private communication). There are 
basically two single-dish and interferometric data combina- 
tion techniques that work well (Heifer et al. I2003> . One is to 
add synthetic short-spacing visibilities in the (m, v) plane (e. g., 
Vogel et al. 1984 . Girart et al. 1999 1. The other is a linear com- 
bination procedure of the single-dish and interferometric maps 
(e.g.. Ye, Turtle & Kennicutt [T99Tl Stanimirovic et al. [T999l . 
We tested these two techniques. 
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Table 2. Parameters of the FCRAO and BIMA Combined Maps 



iri8' -. 









Rms 


Molecular 


Synthesized B 


eam 


Av Noise 


Transition 


HPBW 


P.A. 


(km/s) (JyBeam"') 


HCO+ (l-» 0) 


13'.'3 X ia'3 


-16° 


0.33 0.22 


CS (2^ 1) 


13'.' 1 X ia'3 


-17° 


0.30 0.25 



In order to fill the BIMA visibility hole with synthetic 
FCRAO visibilities, we first deconvolved the FCRAO maps 
from the FCRAO beam pattern using the CLEAN algorithm. 
The beam was modeled as a Gaussian with the FWHM given 
in Table The deconvolved FCRAO map was then multi- 
plied by the interferometer primary-beam response, taken as 
a Gaussian, and centered on the two position observed with 
BIMA (see MGE03). From these two resulting maps, a ran- 
domly distributed sample of visibilities was generated in the 
(m, v) plane of 0-20 light-ns. The number of visibilities gen- 
erated was chosen to have a visibility density similar to the 
BIMA data for the inner part of the (m, v)-plane. To check the 
relative flux calibration of the FCRAO data with the BIMA 
data, FCRAO visibilities were generated in the range of 21 to 
30 light-ns at the same location as the BIMA visibilities. The 
two data sets were compared by vector-averaging the visibili- 
ties in this range over concentric annuli. The values obtained 
from the two instruments agreed to within 15%. 

The linear combination of the maps consists of first gener- 
ating mosaicked BIMA dirty maps, Mbima; as well as the BIMA 
dirty beam, Bbima; then, the "dirty" combined map and beam 
(A^comb, ficomb) are generated as follows (Stanimirovic et al. 
rT9^ : 



comb 



^comb 



I + a 

^bima ^ -^fcrao 

I + a 



(1) 



(2) 



A^fciao is the FCRAO map regridded to match the Mbima cell 
size and map size and weighted with the primary gain of the 
mosaicked BIMA maps; Bfaao is modeled as a Gaussian with 
the FWHM given in Table[n / is the calibration factor, taken as 
1 .0 given the agreement of the relative flux calibration between 
the two instruments; and a is the ratio of the beam areas for 
the FCRAO and BIMA. We adopted a value of a = 0.050 and 
0.042 for the CS and HCO+ data, respectively. The different 
Mcomb were deconvolved from Bcomb using the CLEAN algo- 
rithm. Note that modelling the FCRAO beam as a Gaussian is 
a good approximation: typical sidelobes in the beam pattern of 
a single-dish radio telescope are of the order of -20 dB, which 
will give an eiTor in the combined beam pattern of less than 
lO""', well below the dynamical range of the map. 

As pointed out by Heifer et al. (■2003J . the method of adding 
synthetic single-dish visibilities in the (m, v) plane is parameter- 
dependent, in particular on the number and weighting of the 
synthetic visibilities. Choosing the right combination of pa- 
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ii'is- 
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o 
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iri2' 



19''2l'"00?0 40f0 
a (J2000) 

Fig. 1. Contour maps of the integrated emission over the vlsr 
range 6.05-8.45 km s ' of (from top to bottom) HCO+ 
(1^ 0), H13CO+ (1^ 0), CS (2^ 1), and C^^s (2^ 1). The 
HCO^ contour map is overlaid on the J band image from the 
2MASS All-Sky Survey. The filled triangles show the submm 
sources reported by Visser, Richer & Chandler ^092 (see la- 
bels in the bottom panel). The crosses show the positions of 
infrared sources from the IRAS Serendipitous Survey Catalog. 
From west to east: IRAS 19182-1-1112, IRAS 19183+1109, 
IRAS 19186-Hl 106 and IRAS 19187H-1 107. Contour levels are 
2, 3, 4, . . . 8 times 0.24 K km s 'for CS and HCO+ and 1, 2, 
3, ... 6 times 0.12 and 0.086 K km s ' for C^^S and Hi^CO^, 
respectively. 
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19''2l'"00!0 40'0 
a (J2000) 



Fig. 2. Grey scale image of the optical depth of HCO^ (1^0) 
(top panel) and CS (2— > 1) (bottom panel) derived from the 
integrated maps of Fig. ^ The scale levels are shown on the 
right side of each panel. The overlaid contour maps show the 
integrated emission of HCO^ and CS. The dashed line shows 
the gain at the 0.5 level of the BIMA primary beam response. 
The symbols are the same as those in Fig.[l] 



rameters, we found results almost identical to those of the lin- 
ear method technique. Yet the linear method is much faster 
and simpler, so the combined maps presented here were ob- 
tained using the linear method. Details of the combined maps 
are shown in Tabled 

3. Results and Analysis 

3.1. Spatial distribution of large scale molecular gas 

Figure^shows the integrated emission maps of the four species 
observed with the FCRAO. The maps show that the high- 
density molecular emission arises from a filamentary struc- 
ture oriented in the NW-SE direction. Despite sharing the same 
filamentary structure, the four maps show clear morphologi- 
cal differences. First, the CS is more spread out than HCO^: 
the FWHM along the axis perpendicular to the filament major 
axis is about 2'-3' for HCO^, whereas is about 3'-5' for CS. 
Second, the HCO^ and H'-^CO^ relative intensity peaks do not 
coincide. In particular, at the center of the map they are anti- 



Table 3. Selected positions 



lU. 


Position 


llltip L-UUllLCipal I 


a (J2000) 


6 (J2000) 


a 


19''20"34.3'^ 


11°18'05.5" 


IRAS 19182+1112 


b 


19'"20™43.8'^ 


11°16'05.5" 




c 


19''20™47.9'^ 


11° 16' 05. 5" 


SMM 6 


d 


19''20™54.4'^ 


11°15'25.5" 


N 


e 


19''20"54.7'^ 


11°14'45.5" 




f 


19''20"47.9'^ 


11°14'05.5" 


SMM 3 


g 


19''20"52.0'^ 


11°13'45.5" 


SMM 5 / S 


h 


19''21"'02.9'^ 


11°11'45.5" 


SMM 8 



" Positions N and S correspond to positions selected for study in 
MGE03. 



correlated, i. e., the H CO^ maximum is located in a relative 
minimum of the HCO+. Indeed the H'3CO+ and HCO+ inten- 
sities are comparable at this position. This suggests a very high 
optical depth. In order to estimate quantitatively the effect of 
the opacity on the emission of CS and HCO^ we estimated their 
optical depths from the ratio of the integrated emission maps of 
H'^CO+ and HCO+ and from C^'^S and CS, using the isotope 
abundance ratios shown in Table 0] Figure |2] shows the grey 
scale image of the HCO^ and CS optical depth overlaid with 
the contour map of their integrated emission. This figure shows 
clearly that the HCO^ is highly optically thick, with t ^ 7 for 
most of the emitting region and with a maximum of t = 90 at 
the position of L673 SMM 5. On the other hand, the CS emis- 
sion is moderately optically thick (t ^ 1^). 

The extremely high optical depths of the HCO^ (1— » 0) 
emission imply that this line is not a good tracer of the dis- 
tribution of the HCO^ molecule in L673. Thus, we confirm the 
hypothesis pointed out by MGE03 that the high optical depth of 
HCO^ ( 1 ^ 0) is responsible for the lack of agreement between 
the observed BIMA HCO^ column densities and the theoreti- 
cally predicted column densities. Indeed, the maximum HCO^ 
optical depth measured is located within the region observed 
with BIMA. Strong absorption features and very high optical 
depths of the HCO^ (1^0) emission are also observed in other 
molecular clouds (e. g., Langer et al. 1978 , Girart et al. i2000i) . 

From an inspection to the IRAS catalogues within the re- 
gion shown in Fig. [2 we found four infrared sources. IRAS 
19182-1-1112 and IRAS 19183+1109 are detected only at 
12 and 25 ;um. On the other hand, the other two sources, 
IRAS 19186+1106 and IRAS 19187+1107 are detected only 
at 100 ixm. There are also four weak submillimeter sources 
(Visser, Richer & Chandler I2002> . Only the strongest submm 
source, L673 SMM8, is associated with an infrared source, 
IRAS 19187+1107. All the submm sources ai-e located in 
strong-molecular-Une-emitting regions. Indeed, they are well 
correlated with the relative intensity peaks of H'-'CO^ (1— » 0). 
Visser et al. (2002) classify these submm sources, except 
SMM3, as starless cores. SMM3 is classified as a candidate 
protostar (where there is some evidence of high-velocity CO 
emission). 



Table 4. Physical parameters 
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id. 


Tmb 


Tmb 


V 


^V'' 


Area 




N'' 






(CS) 




(CS) 


(CS) 


(CS) 


(CS) 


(CS) 






(K) 


(K) 


(km s-') 


(km s-') 


(Kkms-') 




(10'^ cm-2) 




a 


1.34+0.14 


< O.W 


+7.0 


0.7 


1.03 + 0.06 


< 1.72 


0.95-1.53^ 




b 


1.64±0.16 


0.32 + 0.06 


+7.2 


0.8 


1.48 + 0.06 


4.21 


3.60 + 0.14 




c 


1.55+0.17 


0.26 + 0.04 


+7.2 


0.8 


1.31+0.05 


3.61 


2.89 + 0.10 




d 


1.45+0.18 


0.24 + 0.06 


+7.4 


0.9 


1.37 + 0.05 


3.53 


2.97 + 0.10 




e 


1.14+0.13 


0.19 + 0.06 


+7.5 


1.1 


1.38 + 0.09 


3.48 


2.97 + 0.19 




f 


J 74+0 21 




+7.2 




1 (i(S + D 0(1 
1 .uo in u.uu 


j.yy 


j.oy X u. 1 J 




S 






+7.2 


u.o 


1 + 0(S 
1 .JO in u.uu 


S 04. 


^81+017 
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4-7 9 


U.J 


n oe -1- n O'l 
y).yo ic u.uj 


9 04. 


1 S(i -1- 07 
i . JO It u.u / 




id. 


Tmb 


Tmb 


v 


^v'' 


Area 






A' 




(HCO+) 


(H'3C0+) 


(H'3C0+) 


(H'3C0+) 


(H'3CO+) 


(H"CO+) 


(H'3C0+) 


(HCO+) 




(K) 


(K) 


(kms-') 


(kms-') 


(Kkms-') 




(10" cm-2) 


(10'3 cm-2) 


a 


1.49 + 0.23 


0.46 + 0.12 


+7.0 


0.4 


0.22 + 0.03 


0.37 


3.86 + 0.59 


2.40 + 0.36 


b 


1.53+0.19 


0.42 + 0.06 


+7.1 


0.5 


0.23 + 0.03 


0.32 


3.91+0.44 


2.42 + 0.27 


c 


0.95 + 0.18 


0.49 + 0.07 


+7.2 


0.7 


0.37 + 0.02 


0.73 


7.16 + 0.47 


4.44 + 0.29 


d 


0.70 + 0.16 


0.34 + 0.05 


+7.4 


0.7 


0.25 + 0.03 


0.66 


4.77 + 0.48 


2.95 + 0.30 


e 


0.54 + 0.15 


0.21 + 0.07 


+7.5 


0.5 


0.12 + 0.02 


0.48 


2.12 + 0.38 


1.31+0.24 


f 


1.17+0.19 


0.78 + 0.04 


+7.3 


0.5 


0.38 + 0.02 


1.10 


8.12 + 0.43 


5.03 ± 0.26 


g 


0.81 +0.17 


0.89 + 0.05 


+7.3 


0.5 


0.45 + 0.02 


> 4'' 


> 18.2 + 0.83 > 


11.3 + 0.51 


h 


2.80 + 0.25 


0.42 + 0.97 


+7.4 


0.6 


0.25 + 0.03 


0.16 


4.10 + 0.48 


2.55 + 0.30 



" The error in the line velocity is lower than 0. 1 km s . 

The error in the line width is lower than 0. 1 km s" ' . 
" Using an abundance ratio CS/C'''S=20 and Tex = 4 K. 

'' Beam averaged column density calculated from the formula in Table 3 of MGE03. 
' Upper limit taken as 3<t, where cr is the sensitivity per channel. 
^ Range of values for t « 1 and t = 1.72. 

« using an abundance ratio HCO+/H"CO+ = 62 (Langer & Penzias fT993l and T^^ = 4 K. 
* We adopt a value of r = 4 as a lower limit. 



By analyzing the kinematical properties of the CS and 
H'^CO^ emission, we found a smooth velocity gradient along 
the major axis of the filamentary molecular structure: from NW 
to SE the v'LSR goes from ~ 7.0 to 7.4 km s"'. In addition, 
east of the center of the map the CS emission is redshifted to 
~ 7.7 kms-'. 

3.1 .1 . Comparison of the emission of the CS (1 0) 
and (2^ 1) lines 

We have also compared the CS (1— > 0) Une measured with the 
Yebes telescope (Morata et al. 19971 with the CS (2-> 1) line 
obtained with the FCRAO telescope, after convolving the spec- 
tra of the (2— > 1 ) line with the 1 .9' beam of the Yebes telescope. 
We compared the resulting spectra at the position of the CS 
(1— > 0) peak, corresponding to the North field of the BIMA 
map, and the position of the peak NH3 intensity, corresponding 
to the South field. The derived Tgx at the North position is the 
same for both lines, 4.2 K. In the South position, where we did 
not observe the C^'*S (1^0) transition, the excitation temper- 
ature derived from the (2— > 1) line, 4.4 K, is very similar to the 
value found in the northern region. 



3.1 .2. Column densities 

As in MGE03, we selected several positions in the region 
mapped with the FCRAO observations in order to obtain the 
physical parameters of the gas and to compare them with the 
chemical model. Table|3lshows the coordinates for each posi- 
tion and their map counterparts, if any. These positions were 
selected mainly to coincide with the 4 positions we selected 
in MGE03, and with the SMM and IRAS sources of the re- 
gion. We also selected two positions not related to any object 
in order to see if they traced "younger" gas. Tablel^shows the 
line parameters for the CS (2— > 1) and H'^CO^ lines, and the 
physical parameters obtained for each position. We adopted a 
value of Tex = 4 K for our calculations, based on the result we 
had obtained from the comparison between the CS (1-^ 0) and 
(2^ 1) hnes (Morata et al. ll997l and MGE03, respectively). 

Tablelsjshows the ratio between the derived column densi- 
ties for CS and HCO^ obtained for each position, which can be 
compared to the results of our chemical model (MGE03). We 
observe that there is a similar trend to the results we obtained 
in MGE03. In the positions where there are no SMM or IRAS 
sources (positions b, d, and e) the CS column density is higher 
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Table 5. Ratio of column densities 



id. MHCO+VMCS) 

a 1.6-2.5° 

b 0.67 ± 0.08 

c 1.54 + 0.11 

d 0.99 ±0.10 

e 0.44 ± 0.08 

f 1.29 ± 0.08 

g > 2.96 ±0.19 

h 1.63 ±0.21 



" The value 1.6 is obtained for t = 1.72. The value 2.5 is obtained 
when T « 1 . 



Table 6. Calculated CS abundances at the positions of the 
SMM sources 



Beam 






N(CS) 




(arcsec) 


(cm-2) 




(cm-2) 


X(CS) 


SMM3 60 


4.3x10^' 


0.15±0.02 3.2 


2.74x10''' 


6.4x10-' 


SMM5 80 


4.0x10-' 


0.18±0.02 3.8 


2.66x10'^ 


6.7xl0~' 


SMM6 100 


2.5x10-' 


0.15±0.02 3.2 


2.56x10'^ 


1.0x10-** 


SMM8 110 


3.8x102' 


0.11±0.02 2.0 


1.08x10'^ 


2.8x10-' 



" From Visser et al. 120(51 



than the HCO^ column density, which would indicate that the 
clumps are chemically young. In the positions where there is a 
SMM or IRAS source, the column density fo HCO^ is higher 
than for CS. This is what we expect for later times, when the 
clumps are more chemically evolved, HCO^ is more abundant, 
and CS starts to suffer freeze-out. We found the highest ratio 
for SMM5 (position g), which we previously classified with the 
BIMA observations as one of the two chemically more evolved 
points, and where the N2H^ emission is more intense. Thus, we 
confirm our previous classification. 

3.1 .3. CS abundances 

Visser et al. ( 2002 1 calculate the masses of the SMM sources 
in L673 from the total flux measured in an area of radius r. 
They also calculated the H2 column density over this area. We 
have convolved the spectra of the CS (2— > 1) and C"''*S (2— > 1) 
emission with a Gaussian with FWHM equal to the r derived 
for each SMM source, in order to estimate the CS column den- 
sity for each SMM position. Table |6l shows the averaged CS 
abundance obtained from the H2 and CS column densities. The 
derived values are in the range 2-10 x lO""^ cm"^. We found 
the lowest value near the IRAS source at the SE part of the 
large scale maps, and the highest abundance at the position of 
SMM6, a submm source that seems to be less centrally peaked 
than other submm sources (Visser et al. I2002> . which would 
indicate that it is in a less evolved stage. 



3.2. Spatial distribution at high angular resolution 

Figurel^shows the composite of the combination process of the 
BIMA and FCRAO observations for the CS (2^ 1) and HCO+ 
(1^0) emission. 

The HCO^ combined map shows a clear filamentary struc- 
ture connecting several clumps, although some of them, par- 
ticularly to the SE, are far from the phase center of the map, 
which makes it difficult to distinguish them from another fila- 
ment. However, as indicated by Fig.|2] the region observed has 
a very high optical depth, so the emission shown in the com- 
bined map cannot be used to describe the true morphological 
and kinematical behaviour of the HCO^ emission. 

On the other hand, the combined map of the CS emission, 
although it has a moderately high opacity (2-4), can be used to 
trace the CS distribution. When we recover all the flux, we find 
that all the most prominent clumps found in the BIMA obser- 
vations are also detected but with more intensity. We believe 
that these clumps, which lie near the phase center of the map, 
are real. We also recover emission from a great deal of other 
structure which was not, or only marginally, detected previ- 
ously. This confirms the predicted effect of the BIMA filtering 
in MGE03. The overall aspect is also of a filamentary struc- 
ture connecting several intense clumps. One of these clumps 
coincides with the nominal position of the source SMM 5. 

The clumps and the filamentary structure of the emission 
are more clearly revealed in Fig. 0] which shows the channel 
maps of the combined BIMA and FCRAO maps for the CS 
(2—> 1) line. Here we find a very complex structure, with sev- 
eral clumps with sizes of ~ 13"-25", appearing in two or three 
of the channel maps. Some of these clumps might coincide with 
clumps found in the HCO^ map. 

3.3. Clump Analysis 

Since the combined high-angular-resolution maps show a high- 
ly clumpy CS spatial and velocity structure, we carried out the 
clump analysis technique developed by Williams, de Geus & 
Blitz (.1994) in our combined BIMA and FCRAO CS (2^ I) 
data cube. This technique consists mainly of two algorithms, 
CLFIND and CLSTATS. CLFIND, which does not assume any 
clump profile, works by contouring the data at a multiple of 
the rms noise of the data cube, then searches for peaks of 
emission to locate the clumps, and then follows them down to 
lower intensities. CLSTATS is a task that calculates integrated 
intensities, peak positions, sizes, velocity dispersion, etc, for 
the clumps found by CLFIND. A very detailed explanation of 
this analysis technique is found in Williams et al. (1994), who 
also found that the CLFIND algorithm obtains the best results 
when the contour interval and the lower contour are set to 2cr. 
However, they recommend using a higher value of the lower 
contour if there is background emission. This is the case for 
the combined L673 CS (2^ 1) maps, where there is some low- 
level (~ 2cr) extended emission. Thus, we run the CLFIND 
and CLSTATS algorithms for the L673 CS channel maps with 
the contour interval set to 0.48 Jy beam"' (2cr) and the lowest 
contour set at 0.96 Jy beam ' (4cr). 
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Fig. 3. Composite of the combination process for CS (2— > 1 ) (bottom panels) and HCO^ (1^0) (top panels). All maps show the 
integrated emission over the v'lsr range 6.05-8.45 km s Left panels: FCRAO maps. Middle-left panels: FCRAO maps cor- 
rected for the primary beam of the mosaicked BIMA maps. Middle-right panels: mosaicked BIMA maps. Right panels: Combined 
BIMA and FCRAO maps. For the FCRAO maps, the contour levels are 3, 4, 5, . . . 12 times 3.17 and 2.80 Jy beam ' km s~' for 
CS (2— > 1) and HCO^ (1^ 0), respectively. For the BIMA and combined maps, the contour levels are 2, 3, 4, . . . 9 times 0.27 
and 0.24 Jy beam ' km s ' for CS (2— > 1) and HCO^ (1^ 0), respectively. The inner and outer dashed lines show the 0.5 and 
0.25 level of the BIMA primary beam response. The symbols are the same as those in Fig. 1. 



A total of 15 clumps were found, all of them with a peak 
intensity of at least ~ Act. Table lists the physical param- 
eters derived for the 15 clumps found by the algorithm. All 
the clumps are clearly resolved at our angular resolution, and 
have deconvolved diameters in the 0.03-0.09 pc range. The ve- 
locity dispersion of the clumps (i. e., the FWHM corrected for 
the spectral resolution) is not very different for the different 
clumps, 0.31-0.65 km s"'. The clumps have masses between 
0.02 and 0.36 Mq. The estimated clump densities range be- 
tween 0.6 and 4.6 x 10"* cm"^. This value of the averaged 
density of the clumps implies that, since the CS (2— > 1) crit- 
ical density is 5 x 10^ cm"^, they are subthermalized. It also 
explains why the derived excitation temperature, 4 K, is signif- 
icantly below the kinetic temperature expected in the cloud, 10 
K. 

Figure |5l shows the location of the clumps derived by the 
algorithm as circles. The areas of the circles are proportional 
to the masses of the clumps. This figure shows a segregation of 
clumps between the northern and southern section of the map. 



The clumps with peak intensities in the 7.2 and 7.5 km s~ 
channels arise from the northern region, whereas the clumps 
that have a peak in the 6.9 km s ' channel are located in the 
southern region. In addition, the southern part contains fewer 
clumps than the northern part, but it is dominated by the largest 
and most massive clump. Interestingly, this clump is associated 
with the two NaH^ clumps found by MGE03. This confirms the 
idea that both N2H^ and NH3 trace the more evolved, massive 
clumps. 

4. Discussion 

4.1. The large-scale structure of L673. Correlation of 
HCO^ and A/2/-/+ emission with dust 

As we have pointed out, the H'^CO^ seems to be well corre- 
lated with the dust as detected by the submillimeter observa- 
tions. This suggests that H'-'CO^ traces the more massive or 
more evolved clumps. At the same time. Fig ^ shows that the 
peak H'^CO^ emission is also coincident with the position of 
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Fig. 4. Channel maps of the CS (2^ 1) emission. The contour levels are -3, 3, 4.5, 6, . . . 12 times 0.24 Jy beam"' (the rms noise 
of the maps). The inner and outer dashed lines show the 0.5 and 0.25 level of the BIMA primary beam response. The symbols 
are the same as in Fig. 1. 



the BIMA peak of N2H+, and a SMM source (SMM 5). In our 
chemical models, HCO^ and N2H+ are late-time molecules. 
Thus, the apparent discrepancy found in MGE03 between the 
column densities observed with BIMA and the expected values 
from the chemical model was due to the extremely high opacity 
of the HCO^ emission. H'^CO^ is found where N2H^ is also 
found, although a little more extended, as can be expected from 
the models, because it is formed a little earlier. Moreover, the 
emission of H'-^CO^, a late-time molecule, coincides with the 
position of the dusty submm sources, which also agrees with 
the theoretical predictions (Taylor et al. il998j . 



4.2. The small-scale structure of L673 

4.2.1 . CS clump distribution properties 

Figure |6l shows the mass and size distribution of the derived 
clumps. It is clear that with only 15 clumps we cannot pur- 
sue a statistical analysis of the mass distribution of the clumps. 
However, with the data available, the mass clump distribution 
is compatible with the low-mass end of the mass distribution 



observed in low-mass star formation clouds, such as in p Oph 



(Motte et al. 119981 Johnstone et al. ■2000J . where the slope of 
dN/dM is not so steep (see Fig.|6j. 

The derived clump masses are significantly below the virial 
mass (see Table 0. This strongly supports the idea that these 
clumps must be transient (Taylor et al. 1996l ll9^ . Recently, 
a detailed study of the illuminated clumps ahead of the HH 2 
object also suggested that clumps in molecular clouds are tran- 
sient (Girart et al. l2002l Viti et al.|2003l|. Fig.|7]also shows that 
the mass of the clumps becomes closer to the virial mass when 
they get bigger and more massive, which is when they have a 
higher chance to form stars. However, the mass calculated for 
the clumps is derived assuming constant CS abundance. If, as it 
seems, the CS abundance is lower for more evolved clumps, the 
slope of the relation shown in Fig. would be steeper, and the 
clumps would probably have a mass closer to the virial mass. 

4.2.2. Largest clump and N2H+ clumps properties 

Figure |8] shows the comparison of two channel maps for the 
combined CS (1— » 2) and N2H+ (1^0) emission. We can 
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Fig. 5. Combined BIMA and FCRAO CS (2^ 1) contour maps for the three channels with the strongest emission (from Fig.|4}. 
The thick circles are centered at the position and at the channel at which each clump shows its highest intensity, as found by the 
CLFIND algorithm of Williams et al. (.1994 ) and described in Sect. 13.31 The areas of the circles are proportional to the masses of 
the clumps. The numbers correspond to the identifier in the first column of Table0 



clearly see that the peak of the CS emission is shifted with re- 
spect to the peak of NaH^. Moreover, the CS emission is more 
spread than N2H^, so we are detecting chemical differentiation 
inside the clump. 

The largest CS clump is not only the clump with the highest 
Mciump/Mviriai ratio (and thus with the highest probability of 
undergoing star formation) but it is also the only one associated 
with NaH^ emission (see Sect. 13. 3> . This confirms the Taylor 
et al. J19961 11998 1 scenario. Interestingly, the density of this 
clump, as derived from the CS emission, is similar to that of 
the other clumps. This may be explained if the emission of the 
N2H^ traces the highest densities of the clump, high enough for 
CS depletion to occur, as has been observed in several dense 
starless cores (e. g., Tafalla et aL EUU^EUn^ . 

We have calculated from the BIMA data the mass traced by 
the N2H^ for this clump. Assuming a N2H+ fractional abun- 
dance of 8 X 10 which was obtained using the H2 column 
density determination from the submm sources of Visser et al. 
J2002> . we could derive a total mass for the biggest clump of 
1.2 Mq. This mass is significantly higher than those of the 
other clumps. In addition, McUmp/Myii- is 0.52, so it is close 
to collapse. 

4.2.3. On the small scale structure of L673 

The previous results, along with the results from the lower- 
resolution analysis in Sect. 13.1.21 show at least the differ- 
ences between the two regions within the BIMA observations. 



T 




Fig. 6. Top Mass distribution of the L673 CS clumps within 
the primary beam of the BIMA maps. The lines show the 
clump mass spectrum for a power-law distribution, oc M 
Bottom Size (equivalent circular radius) distribution of the 
L673 CS clumps. The dotted line shows the radius of the CS 
(2^ 1) synthesized beam. Note that all the clumps are well re- 
solved at our angular resolution. 
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Table 7. Clumps within the BIMA primary beam 
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Optical depth derived at the peak position of the clump from the 
intensity ratio of the CS (2^ 1) and C^^'S (2^ 1) FCRAO maps. 
FWHM of the line corrected for channel resolution as described 
in Appendix A of Williams et al. (1994) 

Equivalent circular radius corrected for beam size as described in 
Appendix A of Williams et al. (1994) 

Clump mass derived from the integrated emission obtained with 
CLSTATS corrected for the CS (2— » 1) optical depth, t, and as- 
suming a constant excitation temperature of r„ = 4 K. An abun- 
dance ratio of X[CS/H2] = 10"** is adopted. 
Virial mass obtained from Myj, = 5 ARcrl/aG with cr,, = 
Av/2.355 and assuming an inverse-square power law density pro- 
file (a=5/3). 

Average volume density obtained assuming a sphere of mass 
and radius AR. 



The northern region shows a gas less chemically evolved, and 
smaller and more numerous clumps. The southern region con- 
tains the more evolved gas and the biggest clump. 

The clumps detected in our observations seem to be of a dif- 
ferent kind than the ones usually studied in surveys of "starless 
cores" (Tafalla et al. 2002; Bergin et al. l2001> . Those are usu- 
ally older cores, which have already achieved the more evolved 
and more massive stage of evolution, as is clearly indicated by 
the detection of late-time molecules, as NiH^ and NH3, and by 
the high level of depletion of some molecules onto dust grains. 
We find one such clump in our region, which might be just 
arriving at that evolved phase. This clump is at the same time 
the most massive, shows emission of N2H^, and different emis- 
sion peaks for CS and N2H^, which could point to depletion. 
However, we detect many other clumps, less massive, that show 
almost no emission of late-time molecules, and which probably 
will disperse before they can achieve the "starless core" stage. 

We face the question of what the origin of these transient 
clumps could be. They are probably caused by interstellar tur- 
bulence (see e.g. Padoan 1995 ; Ballesteros-Paredes et al."1999 l. 
In support of this idea, we observe that the derived masses 
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Fig. 7. Relation between the clump to virial mass ratio with re- 
spect to the clump mass (top) and the equivalent circular radius 
{bottom). The open circle represents the most massive clump 
when the mass derived from the NaH^ emission is considered. 

for our clumps are around 0.1 Mq, which is also the mass- 
scale of turbulent fragmentation in the ISM (Padoan Il995> . 
This small scale turbulence is probably related to large-scale 
flows in the diff'use Galactic interstellar medium, which sev- 
eral studies (Hartmann et al. 1200 II Bergin et al. l2004> propose 
as a mechanism that can lead to the fast formation and dis- 
persal of molecular clouds, which could explain the apparent 
rapid star formation and short cloud lifetimes observed in the 
regions of the solar neighborhood. On a smaller scale, several 
authors postulate the possible existence of transient clumps. 



Falle & Hartquist J2002> propose that slow-mode MHD waves 
could create, under suitable conditions, large density enhance- 
ments, of the order of a factor of 30, which could last on the 
order of 1 Myr, before dispersing back to the original density. 
At the same time, Vazquez-Semadeni et al. ( 120051 12003> and 
Ballesteros-Paredes et al. (2003,) also argue that cores within 
molecular clouds may not be in hydrostatical equilibrium, and 
would be transient features generated by the dynamical flow 
in the cloud, not necessarily requiring strong magnetic fields. 
These cores must then either collapse, if they have an excess of 
gravitational energy, or re-expand and merge back into the sur- 
rounding molecular cloud. The re-expansion time is expected 
to be longer than the compression time because of the retard- 
ing action of self-gravity, and is estimated to be of the order of 
a few free-fall times. This would be consistent with the result 
that typically there are more starless than star-forming cores in 
molecular clouds and that, as in our case, most of the cores do 
not appear to be gravitationally bound. Only sufficiently mas- 
sive cores would eventually undergo local collapse. 

In any case, such models would have consequences for the 
physical and chemical evolution of these regions (Williams & 
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Viti 12002 1). Firstly, the chemistry would be mainly "young", 
as the molecules would be destroyed when the gas returns to 
a more diffuse state. Secondly, this mechanism would halt the 
process of freeze-out of molecules onto dust, as the frozen-out 
molecules would also return to the gas phase as the density of 
the gas goes back to the diffuse state. Thus, these regions would 
not be completely depleted of molecules in the gas phase, as is 
apparently shown by some observations (Gibb & Little 1998 1. 
Finally, the mass distribution of these transient clumps proba- 
bly determines the low-mass star formation rate. 

However, an alternative possibility is that the more massive 
clumps were formed by aggregation of smaller clumps into a 
massive one, which could then be closer to collapse. The analy- 
sis of our clumps allows us to make some rough estimates to see 
if this can happen since transient clumps are expected to sur- 
vive for ~ 1 Myr (Falle & Harquist 2002; Vazquez-Semadeni 
et al. 2004). Two limiting cases can be studied. One is that 
of a 1 -dimensional structure, that is, the clumps are assumed 
to move inside a filamentary structure, with relative velocities 
(between clumps) of 0.3 km s This is the typical difference 
in the radial velocities of the clumps, which should roughly 
give a reasonable lower limit for the true relative velocities. 
In the course of ~ 1 Myr any of the observed clumps would 
have travelled ~ 63000 AU or ~ 3'.5, which is the size of 
the field of view of our maps. The other case is that of a 3- 
dimensional (3-D) structure with similar scales in the three di- 
mensions. In this case the projected area of this volume in any 
particular direction is Aproj - Vtof , where Aproj and Vtot are 
the area and volume, respectively. Let us assume that in a vol- 
ume Vtot there are n identical spherical clumps of radius R dis- 
tributed randomly, with a random velocity distribution with an 
average velocity {v3d). The 3-D filling factor of the clumps for 
this volume is /jd - nVa/Vtot - "^ciMp^oj' where Vd is the 
clump volume. The volume swept up by a clump in a time t is 
VswiO - v^D tJiR^. For a given clump, the time scale needed 
to cross another clump, tx, is such that /3D Vsw(fx) = Vd- 
From the previous equations, the average crossing scale time 
can be written as (fx) - A^J^^I{n{vT,^) nR^). Since what we 
measure is the radial velocity, Vobs and (v3d) = 2(vobs), then 
{tx) = Ap|!^j/(2n(vobs)7r^^)- For Aproj equal to our field of 
view (~ 3.5 X 2.'0), which includes 15 clumps, and using a ra- 
dius equal to the average value for the measured clumps (from 
Table 7), R - 18", we found that {tx) -6x10^ yr. Then it could 
be possible that star formation was the result not of the gravita- 
tional collapse of just one clump, but of a "merging" effect. In 
this case, the geometry of dense cores would be far from spher- 
ical symmetry, even the density distribution would not follow 
n(r) oc r^^. Tafalla et al. ('2004*1 find that their starless cores are 
not spherical, but that clear inhomogeneities are found in them. 

5. Conclusions 

We have made a multitransitional study of the spatial distri- 
bution of the molecular emission with the 14 m FCRAO tele- 
scope in the starless core found in the L673 region, which had 
been previously observed at lower angular resolution (Morata 
et al.[T997^ and with the BIMA array (MGE03). The main goal 
of these observations was to combine them with the array ob- 




HA (J2000) EA (J2000) 

Fig. 8. Comparison of the channel maps of the combined CS 
(2— > 1) emission (solid contours) for the 7.16 {left) and 6.87 
km s ' {right) channels; and the BIMA N2H^ (1— > 0) emission 
(contours and grey scale) for the -0.8 {left) and -1.1 km s"' 
{right) channels. The contour levels for the CS maps are the 
same as in Fig.|3] The contour levels for the N2H^ maps are 3, 
3.5, 4, 5, 6, and 7 times 0.20 Jy beam"' (the rms noise of the 
maps). The symbols are the same as those in Fig.^ 

servations in order to obtain high resolution maps of the full 
emission, and thus be able to confirm that the clumpy structure 
detected with the BIMA telescope was real. The main results 
were: 

1. With the FCRAO we detected emission in the CS (2-> 1), 
C^'^S (2^ 1), HCO+ (1^ 0), and H'^CO^ (1^ 0) lines. 
Although every map shows clear morphological differ- 
ences, we found that the high density molecular emission 
in the four molecules arises from a filamentary structure 
oriented in the NW-SE direction. The ratio between HCO^ 
and H'^CO^ also confirms that the apparent lack of HCO^ 
(1— > 0) emission in the BIMA maps was due to the ex- 
tremely high optical depth in this line. We also found that 
several submm sources detected in this region are located in 
strong molecular line emitting regions, and that they seem 
to be well correlated with the relative intensity peaks of 
Hi^CO^. 

2. We selected several positions to study the local physical 
parameters. We confirm the results of MGE03: the column 
density of CS is higher than that of HCO^ for positions 
were the gas is chemically younger, i.e. not associated with 
IRAS or submm sources. In the positions where there is 
a submm or IRAS source, the HCO^ column density is 
higher than that of CS. 

3. The BIMA and FCRAO combined maps of the CS emis- 
sion show that all the prominent clumps found in the BIMA 
observations are also detected but with more intensity. We 
believe that these clumps are real. We also recover emis- 
sion from other structure which was previously undetected 
or only marginally detected. The overall aspect is also of a 
filamentary structure connecting several intense clumps. 
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4. We found a total of 15 clumps in our combined data cube. 
All of them are resolved at our angular resolution, with di- 
ameters in the 0.03-0.09 pc range. Their estimated masses 
range between 0.02 and 0.36 Mq, and their densities be- 
tween 0.6 and 4.6x10"^ cm"^, which indicates that the CS 
(2^ 1) emission is clearly subthermalized. 

5. There is a clear segregation in the properties of the de- 
tected clumps between the northern and southern regions 
of the map. The northern region shows the less chemically 
evolved gas, and less massive but more numerous clumps. 
On the other hand, the southern region contains the more 
evolved gas and the most massive clump. This clump is 
also the only one associated with BIMA N2H^ emission. 
Moreover the H^^CO"^ emission correlates well with the 
dust emission traced by the submillimeter observations, and 
its emission peak also coincides with the BIMA N2H+ peak 
and the SMM 5 source. AH these results support the predic- 
tions of our chemical model, which proposed that late-time 
molecules, such as N2H+ and HCO^, should be found in 
the same regions. 

6. The clump masses derived are below the virial mass, which 
means that these clumps must be transient, as only the more 
massive ones are able to condense into stars. 

7. The properties of the clumps we detect seem to be differ- 
ent from those of the "starless cores" recently found in the 
literature. Although we might have detected one such core, 
the rest are in an early evolutionary stage, after which they 
will most likely disperse. 

We find that the simple chemical model proposed in 
MGE03 to explain the properties of the CS core of L673 is 
confirmed by the combination of single-dish and array observa- 
tions. However, there is still a lot of work to be done on the for- 
mation mechanisms of these clumps and on the possibility that 
the merging of smaller clumps can lead to more massive and 
gravitationally unstable clumps. Further observing and mod- 
elling of the chemistry of these clumps and of the interclump 
gas will be needed in order to test the transient nature of the 
clumps. 
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